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We construct a new isentropic equation of state (EOS) at finite temperature “CRover” 
on the basis of the hadron-quark crossover at high density. By using the new EOS, 
we study the structure of hot neutron stars at birth with the typical lepton fraction 
(Y/ = 0.3 — 0.4) and the typical entropy per baryon (5 = 1 — 2). Due to the gradual 
appearance of quark degrees of freedom at high density, the temperature T and the 
baryon density p at the center of the hot neutron stars with the hadron-quark crossover 
are found to be smaller than those without the crossover by a factor of 2 or more. Typical 
energy release due to the contraction of a hot neutron star to a cold neutron star with 
the mass M = 1.4Mq is shown to be about O.O4M0 with the spin-up rate about 14 %. 


Subject Index Neutron stars, Nuclear matter aspects in nuclear astrophysics, Hadrons and 
quarks in nuclear matter, Quark matter 


Introduction: In the core-collapsed Type-II supernova explosion, the proto-neutron star 
(PNS) with the radius ~ 100—200 km is formed. During the first few seconds after the 
core bounce, the PNS undergoes a rapid contraction and evolves into either a “hot” neutron 
star (NS) with the radius ~ 10—20 km or a black hole. The hot NS at birth in quasi¬ 
hydrostatic equilibrium is composed of the supernova matter with the typical lepton fraction, 
Yi = Ye + Yi, ^ 0.3 — 0.4, and the typical entropy per baryon^, 5 ~ 1 — 2: They are caused 
by the neutrino trapping at the baryon density p exceeding lO^^g/cm^. With this as an initial 
condition, the hot NS contracts gradually by the neutrino diffusion with the time scale of 
several tens of seconds and evolves to a nearly “cold” NS with Tj, ~ 0 and S' ~ 0), unless 
another collapse to a black hole does not take place [1-3]. 

The hot neutron star provides us with various information on the properties and dynamics 
of high density matter [4, 5]. The purpose of this Letter is to study the hot neutron star at 
birth with degenerate neutrinos on the basis of an equation of state (EOS) "CRover" at finite 
temperature T which is newly developed on the basis of the hadron-quark crossover picture at 
high baryon densities. Such an EOS for “cold” neutron-star matter has been previously stud¬ 
ied by the present authors [6, 7]: It was shown that a smooth crossover from the hadronic 
matter to the strongly-interacting quark matter around p ~ 3/9o (po = 0.17/fm^ being the 
normal nuclear matter density) can support the cold neutron star with the maximum mass 
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greater than the observed masses of massive neutron stars, M = (1.97 ± 0.04)Mq [8] and 
M = (2.01 ± 0.04)Mq [9]. Similar conclusions for cold neutron stars have been recently 
reported by several other groups [ 10 - 12 ]. 

In this Letter, we will generalize the idea of the hadron-quark crossover to the system at 
finite T by considering the Helmholtz free energy F{N, V, T) as a basic quantity to interpolate 
the hadronic matter and the quark matter. Here N and V are the total baryon number and 
the total volume, respectively. By using our EOS for supernova matter with fixed Yi and S, 
we study not only the bulk observables such as M-R relation of the hot neutron star at birth, 
but also the temperature, density and sound velocity profiles inside the star. We note here 
that the hot neutron stars with the hadron-quark mixed phase has been studied previously, 
e.g. [13-16]. Such a mixed phase generally leads to soft EOS, so that it is rather difficult 
to sustain 2 M 0 NSs. On the contrary, our hadron-quark crossover approach does not suffer 
from the problem, since it leads to stiff EOS in the crossover region. 

EOS for supernova matter at finite T with hadrons and leptons: A set of hnite- 
temperature Hartree-Eock (HE) equations is solved for an isothermal matter composed of 
n,p, e“, 6 "*“, t'e and Pg, under the conditions of charge neutrality, chemical equilibrium and 
baryon and lepton number conservations. It gives the single-particle energies, chemical poten¬ 
tials and mixing ratios of nucleons and leptons, so that the hadron+lepton EOS (HL-EOS) 
for a given set of T, p and Y) is obtained. Then, we calculate the thermodynamic quantities 
such as the internal energy E, the entropy S and the free energy F(= E — TS) (details are 
shown in [17]). We neglect hyperons (Y) in the hadronic matter in this Letter, since they do 
not appear as hadronic degrees of freedom once the hadron-quark crossover takes place at 
around p ~ 3pQ as shown in [ 6 , 7]. 

We employ V)vjv(= Vrsc + t^TNl) for the effective NN interaction in the HE calcula¬ 
tion: Here Vrsc is a two-nucleon potential based on the G-matrix calculations with the 
Reid-soft-core potential: It depends on the total baryon density and the density difference 
between neutrons and protons. On the other hand, Ltni is based on the phenomenological 
three-nucleon potential of Friedman-Pandharipande type [19] and is written in a form of 
p-dependent two-body interaction. We introduce C/tni to assure the correct saturation of 
symmetric nuclear matter with the nuclear incompressibility k = 250 MeV consistent with 
experiments. This is called the TNI2-EOS [17]. The maximum mass of the NS at T = 0 with 
TNI2-EOS alone is 1.63Mq. However, the hadron-quark crossover inside the NSs can bring 
this well above 2Mq as shown in [ 6 , 7]. 

EOS for supernova matter at finite T with quarks and leptons: To obtain the 
quark+lepton EOS (QL-EOS), we follow [ 6 , 7] and apply the ( 2 -|-l)-flavor Nambu—Jona- 
Lasinio (NJL) model, 

r ® 

-Cnjl = - m)q + ^[(gA“g )2 -h (gfysA^g)^] - [detg(l 75)9 -h h.c.] - ( 1 ) 

a=0 

Here the quarks qi [i = u, d, s) with current quark masses rrii have three colors and three 
flavors. The term proportional to in Eq.(l) is a 1/(3)^ x U{3)ji symmetric four-fermi inter¬ 
action where A“ are the Gell-Mann matrices with A'^ = \/2j3 1. The term proportional to G^ 
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is the Kobayashi-Maskawa-’t Hooft (KMT) six-Fermi interaction which breaks U{1)a sym¬ 
metry explicitly. The last term with the positive coupling represents a universal repulsion 
among different flavors. A standard parameter set obtained from the hadron phenomenol¬ 
ogy in the vacuum [18] reads A (UV cutoff)=631.4 MeV, GgA^=3.67, G^A®=9.29, mu^d—5.5 
MeV, and =135.7 MeV. As for the universal vector coupling, we consider the typical 
value, g^jGg = 0.5. The maximum mass of cold NS with hadron-quark crossover with these 
parameters reads 2.6Mq (see Fig. 19 and Table 7 of [7]). 

The Gibbs free-energy 11 is a sum of the quark contribution and the lepton contribution, 
n(^, V, T]fii)= nquark(li, V, T] fii) + niepton(li, V, T ; /i/) with fj. and m being the baryon chem¬ 
ical potential and the lepton chemical potential respectively. ^ The former in the mean-field 
approximation is 


^ quark 


V 



where rii = (gjgi) is the quark number density in each flavor (i = u,d,s). Si is the quark 
propagator, which can be written as 



(3) 


j 


^ Thermal contribution from massless gluons above the crossover density is less than 4% of the 
thermal part of quarks for T < 30 MeV, so that we do not consider gluons in this Letter. 
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Fig. 1 Schematic illustrations of internal structure and relevant EOS in hot and cold NSs with M = 2Mq 
obtained by a new EOS "CRover" with hadron-quark crossover at finite temperature. We take {Yi,S) = (0.3,1) for 
the hot NS in this figure. BPS, HS, TNI2 and NJL are the EOS adopted (details are in the text). 
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where ui = (2£ + l)7rT is the Matsubara frequency and is the effective chemical potential. 
The lepton part Oiepton(M) corresponds to a relativistic and degenerate system of 

electrons, neutrinos and their anti-particles. The charge neutrality and chemical equilibrium 
among quarks and leptons are imposed. The Helmholtz free-energy F{N, V,T\Yi) is obtained 
from the Gibbs free-energy fii) by the Legendre transformation. 

Hadron-Quark crossover at finite T: 

Let us now generalize the idea of the hadron-quark crossover introduced in our previous 
works [6, 7] to the system at finite T. We consider the following interpolation of the Helmholtz 
free-energy per baryon F between the hadron+lepton free-energy per baryon (T^hl) ^ind the 
quark+lepton free-energy per baryon (Tql), 

F(p,r;yz) =THL(p,r;yzV_(p,r)+FQL(p,r;yz)ru+(p,T). (4) 

Here W- and = 1 — W- are the weight functions. The typical temperature we consider is 
30 MeV or less. This is considerably smaller than the thermal dissociation temperature (~ 
200 MeV) of hadrons into quarks and gluons [20], so that we assume W- being T-independent 
and takes the same form as given in [6, 7]: 

w±{p,T) ^ w±{p) = , (5) 

where p and T are the phenomenological interpolation parameters which characterize the 
crossover density and the width of the crossover window, respectively. We take (p, T) = 
{3po,po) as typical values to account for the 2Mq neutron stars at T = 0 [6, 7]. Note that 
hyperons do not appear in the hadronic phase with such a relatively low crossover density. 

From Eq.(4) with Eq.(5), other bulk quantities, the entropy per baryon (S) and the energy 
per baryon (£’), are derived by using thermodynamic relations, S = —dF/dT and E = F + 




Fig. 2 (a) The entropy per baryon 5 as a function of the baryon density p for Yi = 0.3 with T=10 MeV (red), 

15 MeV (blue) and 20 MeV (green). Solid lines and dashed lines correspond to the cases with crossover and without 
crossover, respectively, (b) The temperature T of the isentropic matter as a function of the baryon density p for Yi = 0.3 
with S = 1 (red) and S = 2 (blue). Solid lines and the dashed line correspond to the cases with crossover and without 
crossover, respectively. 
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(C) 

Fig. 3 (a) The isentropic pressure P as a function of baryon density p for (Y^, S) =(0.3, 1), (0.3, 2) and (0.4, 1). 

The black line corresponds to the EOS for cold neutron star matter. The crossover window is shown by the shaded 
area on the horizontal axis, (b) The energy per baryon E with the same set of Yi and S as (a), (c) P as a function of 

£. 


TS, as 


S{p,T;Yi) = Sni.{p,T;Yi)w.{p) + S^i^{p,T;Yi)w+{p), (6) 

E{p,T;Yi) = EnL{p,T-,Yi)w.ip) + EQUP,T-,Yi)w+ip). (7) 

The thermodynamic quantities for isothermal matter with p, T and Yi can be converted 
into those for isentropic matter by means of the T-p relationship constrained by the constant 
entropy per baryon, S =const. Then the isentropic pressure is obtained as 


P{p,T{p)-YpS) 


dE 

dV 



S,N 


dE 

dp 


S,N 


( 8 ) 
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(a) 



Fig. 4 (a) The neutron star mass M as a function of the total baryon number • Red (blue) curves correspond 

to the hot (cold) neutron stars with crossover. The dotted lines correspond to the case without crossover, (b) The 
energy release AE as a function of the cold neutron star mass Mcoid- S) = (0.3,1) is adopted. 


EOS for supernova matter below pq : 

Below the normal nuclear matter density po, we use the thermal EOS which consists of an 
ensemble of nuclei and interacting nucleons in nuclear statistical equilibrium by Hempel and 
Schaffner-Bielich [21]. We call this as HS EOS in the following. Various EOSs below pQ do not 
show quantitative difference as discussed in [22]. Once the baryon density becomes smaller 
than the neutron drip density 10“^/90) temperature becomes smaller than 0.1 MeV, so 
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that we switch to the standard BPS EOS [23]. In Fig.l, the internal structures of hot and 
cold NSs with our isentropic EOS are illustrated. 

EOS for supernova matter with hadron-quark crossover: To carry out the conver¬ 
sion from the isothermal matter with hxed T and Yi to the isentropic matter with hxed S 
and Yi, we consider S{p,T,Yi) as a function of p/po- In Fig.2 (a), we show S for T—10, 15, 20 
MeV and Yi = 0.3 with crossover (solid lines) and without crossover (dashed lines). One hnds 
that the entropy per baryon turns out to increase at high densities, once the quark degrees 
of freedom start to appear. Fig.2 (b) shows the temperature as a function of p/po under 
isentropic conditions for S = 1,2 and Yi = 0.3 with crossover (solid lines) and for S = 2 and 
Yi = 0.3 without crossover (the dashed line). Due to the appearance of the colored degrees 
of freedom, T is suppressed above the crossover density. 

In Fig.3(a), (b) and (c), we show the isentropic pressure P{p,T{p),Yi) and the energy per 
baryon E{p,T{p),Yi), and the P — e relation, respectively, for the three characteristic sets, 
{Yi,S) = (0.3,1), (0.3,2) and (0.4,1). Here e = pE is the energy density. For comparison, 
the EOS of cold neutron star matter (T = 0 without neutrino degeneracy) is also shown by 
the black solid lines. We call this new EOS with hadron-quark crossover as “CRover” and 
will post the numerical tables for different combinations of S and Yi in the EOS database 
(EOSDB) [24]. 

Structure of hot neutron stars: We now solve the Tolman-Oppenheimer-Volkov (TOV) 
equation to obtain the structure of hot neutron star at birth; 

dP C 

— = (M(r)+47rPr3)(e + P)(l-2GM(r)/r)“\ 

pr 

M{r) = / A7Tr'‘^e{r')dr'. (9) 

Jo 

Here we have assumed spherically symmetric stars with r being the radial distance from 
the center. In the following, we will consider the hot neutron stars with typical values, 
(y;,5) = (0.3,1). 

In Fig. 4 (a), we show the relation between the gravitational mass M and the baryon number 
Nb for hot and cold neutron stars with and without crossover. First of all, the crossover leads 
to heavier neutron stars for given Nb due to the stiffening of the EOS. Furthermore, hot 
neutron stars have larger mass than the cold ones for given Nb- Typical energy release for 
■^cold = 1.4M0 due to the contraction reads AE = M^^ot — Tfcoid ~ O.O4M0. In Fig. 4 (b), 
we show AE in the unit of Mq as a function of the mass of cold neutron star, Mcoid- 

In Fig. 5 (a), the mass-radius (M-R) relations of the hot and cold neutron stars with 
hadron-quark crossover are shown. Although the maximum masses are similar between hot 
and cold neutron stars, the stiffening of the EOS at finite T at low density (see Fig. 3 (c)) 
makes the hot neutron star bigger in size, especially for neutron stars with small M. In Fig.5 


^ In the approach based on the hadron-quark mixed phase, the softening of the EOS due to the 
mixed phase is tamed by the finite temperature effect, so that the maximum value of Nb for hot NS 
becomes larger than that of the cold NS. This is why there are hot NSs which eventually collapse 
into black holes [13, 15]. Similar observation in the case of the pion condensation has been reported 
ealier by one of the present authors [4]. On the other hand, in our hadron-quark crossover picture, 
the hot EOS and cold EOS are similar except for the low density region, so that the delayed collapse 
to the black hole does not take place. 
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Fig. 5 (a) Mass-radius relationship for (Yj, S) = (0.3,1). Red: hot neutron stars with crossover. Blue: cold neutron 

stars with crossover, (b) The sound velocity squared uj as a function of the distance from the center of 1.4Mq neutron 
star. Colors on each line are the same as in (a). 


(b), we plot the local sound velocity squared v^(r) as a function of r for M = 1.4M0 obtained 

by 


2/ V 

Vsip-,yi,s) = 


Yi,S 


dP{p,T{p)-Yi,S)/dp 

de{p,T{py,Yi,S)/dp 


( 10 ) 


with p{r) obtained by the TOV equation. The sound velocity for cold neutron star is larger 
(smaller) at higher (lower) density than that of the hot neutron star. 

In Fig. 6(a) and (b), the central temperature Tcent and the central density Pcent- of the hot 
neutron stars with and without the crossover are plotted as a function of M. The temperature 
decreases due to the appearance of the quark degrees of freedom (see Fig.2(b)), so that the 
central temperature becomes significantly smaller with crossover as shown in Fig. 6(a). The 
central density of the star becomes significantly smaller with crossover as shown Fig. 6(b) 
due to the fact that the EOS is stiffer with crossover. 

To compare the internal structure of the hot neutron star with and without the crossover in 
terms of the temperature and density, we plot T(r) and p{r) in Fig. 7 (a) and (b), respectively. 
Here we consider the hot neutron star with a canonical mass M = 1.4M0. Under the presence 
of the hadron-quark crossover, the EOS is stiffer and hence the central density is smaller, so 
that the temperature and density profiles of the star are more uniform as compared to those 
without the crossover. 


Summary and concluding remarks: In this Letter, we have discussed the properties 
of hot neutron stars at birth on the basis of a new EOS “CRover” for supernova matter 
with hadron-quark crossover. Such a crossover leads to the EOS stiff enough to sustain 2Mq 
neutron stars. A noticeable point is that the crossover plays important roles not only to 
generate the stiff EOS but also to lower the internal temperature of hot neutron stars. Such 
suppression of temperature originates from a combined effect of the isentropy nature of the 
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Fig. 6 (a) The central temperature Tcent as a function of M. (b) The central density pcent as a function of the 

neutron star mass M. The solid (dashed) lines correspond to the EOS with (without) crossover, and {Yi,S) = (0.3,1) 
is taken. 



(a) (b) 


Fig. 7 (a) The temperature profiles of the hot neutron star with M = IAMq and (Y^, S) = (0.3,1). Solid (dashed) 

lines correspond to the EOS with crossover and without crossover respectively, (b) The density profiles of the same 
neutron star as the case (a). 


supernova matter and larger entropy for given temperature due to the quark degrees of 
freedom. Given baryon number, hot neutron stars have larger radius and larger gravitational 
mass caused by the high lepton fraction and the thermal effect. This suggests that, during 
the contraction from hot to cold stars, gravitational energy is released and simultaneously the 
spin-up takes place. In the present study, the released energy is about 0.04 Mq and the spin- 
up rate is about 14 % (assuming the conservation of angular momentum) for M^oid = 1-4:Mq 
of evolved cold neutron stars. It is worth noting here that the maximum baryon number 
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-^max is smaller for the hot neutron star if there is a crossover. This implies that there could 
by massive neutron stars which can be reached only by the mass accretion after the birth. In 
the present study, we have used temperature-independent interpolation function to construct 
the EOS with crossover. Although this is justified for the supernova matter with temperature 
less than 30 MeV, further investigations would be necessary to have the EOS applicable to 
wider range of temperature in supernovae and neutron star mergers. 
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